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analysis is implemented by introducing a fault in the network. Here, a power 
IPSL , fault is injected at bus 4 between the time interval of 2 to 3 sec to analyse the 
Load flow analysis stability of the system. The voltage stability of the system is analysed for the 
Power system stabilizer network with and without automatic voltage regulator (AVR). The AVR unit 
is tuned by using power system stabilizer (PSS). The results are examined by 
simulating the network using openmodelica connection editor. From the 
simulation results the per unit voltages and angles at all 10 buses are 
determined for the network with and without AVR. By comparing both the 
results it is proved that the network with AVR has better voltage stability 
than the other. Thus, the voltage stability of the system is improved by 
connecting the generator with AVR and PSS. 
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1. INTRODUCTION 

In this modern civilization, the electrical utilities need to be automated to fulfill the requirements of 
human activities. For automation various control processes are implemented to control the utilities. These 
processes may be of linear or non-linear in nature which results in reducing the quality of power. In current 
situation the demand for electricity increases like the increase in demand for gold in our country. So, we 
should ensure that all the end users are supplied with quality power [1]-[3]. One can be confident about the 
power quality when there is enough knowledge about stability of voltages, different bus angles and power 
flows in a network. Hence, it is necessary to analyze power or load flow studies of network (either 
transmission or distribution) to enhance its efficiency and quality [4]-[6]. Depends on the power flow studies 
the following buses are framed, a) slack or reference bus, b) generator bus, and c) load bus 

Here, the slack bus is provided with both voltage magnitude and angle. So, it is also called as 
reference bus in the network. The generator bus is called as PV bus because the real power and voltage 
magnitude are specified for it [7], [8]. The load buses they are generally known as PQ buses [9], [10]. The 
per unit magnitude and phase angle of voltage for load buses are determined by using numerical methods 
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[11]. The numerical methods include Newton’s Raphson method, Gauss Seidel Method, Fast decoupled 
method etc [12], [13]. In this research work, a well-conditioned power flow problem is framed and solved 
using Newton-Raphson method. In NR technique, the non-linear equations are solved by polar coordinate 
method. [14]-[16]. 

In general, NR method follows the steps, 

Step 1 : An admittance matrix is framed. 

Step 2 : Initial voltage magnitude and phase angles are assigned to all the load buses connected in the 
network. 

Step 3: Calculate P and Q for all load buses. 

Step 4 : Frame a Jacobian matrix. 

Step 5 : Determine AVi and Aoi by using Jacobian matrix. 

Step 6 : Update the voltage magnitude and phase angle for all the load buses. 

Step 7 : Continue with consecutive iterations until the real and reactive powers are within the specified 
range [17]-[19]. 

The open modelica software is an open-source software and the circuit design is purely dependant 
on the components from Itesla Power System Library. Since, most of the loads used today are of non linear in 
nature, they inject harmonics in the distribution network [20]-[22]. These harmonics results in voltage 
disturbances in the system. Though several methodologies like implementing custom FACT devices, filtering 
techniques are introduced to reduce the harmonics, load flow and stability analysis are the prime significant 
methodologies to perform detailed analysis of per unit voltages at all buses [23]-[25]. 

The paper is organized in such a way that Section 2 explains the proposed 10 bus loop distribution 
system. Section 3 describes the simulation results of load flow analysis of the distribution network with and 
without AVR and PSS. Also, comparison of the simulation results is tabulated in Section 4. Section 5 
concludes the power flow analysed for the networks. 


2. PROPOSED 10 BUS LOOP DISTRIBUTION NETWORKS 

The single line diagram of 10 bus loop distribution network is given in Figurem1. This network is 
energized by using a single generator. It is a 100 MVA base system. In this diagram Bus 0 is a PV bus which 
is of 11kV with active power as Ipu and voltage magnitude as Ipu. Buses from 1-10 are load buses. The 
active power rating is Ipu and reactive power rating is 0.5pu. Table 1 describes the power line parameters of 
the distribution network. By using the above parameters, the voltage magnitude and angles of load buses are 
calculated using NR method. The steps and equations involved in analysis is given: 
Step 1 : Initially, the voltage magnitude of load buses is assumed as 1 and phase angle as 0. 
Step 2 : Calculate P and Q using the (1) to (4). 


[p.ir; p.11] = -[sin(delta), cos(delta); -cos(delta), sin(delta)| *[id; iq] * I (1) 
[p.vr; p.vi] = [sin(delta), cos(delta); -cos(delta), sin(delta)|] *[vd; vq] * V (2) 
P = -p. vr * p. ir - p. vi * p. i (3) 
Q = -p.vi * p. ir + p. vr * p. ii (4) 


Step 3 : Jacobian matrix is formed. 
Step 4 : Change in voltage magnitude and angle of load bus is determined using (5) and (6). 


AV = sqrt (p. vr^2 + p. vi*2) (5) 

Ad = atan2 (p.vi, p. vr) (6) 
Step 6 : The updated value of Voltage and angle is computed by (7) and (8). 

Vnew = Vold + AV (7) 


onew = dold + Ad (8) 
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Step 7 
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: The process is repeated from step 2 for next iteration. The iteration is continued until the load 


buses have minimal error, which is within the specified tolerance level. After completing sufficient iteration, 
magnitude and angle of voltage in all PQ buses are concluded. 


3. 
3.1. 


Power Fault at 
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Figure 1. Single line diagram of 10 bus distribution network 


Table 1. Power line parameters 


S.NO LINEDATA VALUE 

1 Inductance of line 1.98mH 

2 Resistance of line 0.369 

3 Capacitance of line 25Mf 

4 Length of the Distribution Feeder 10km 

5 Line Impedance 0.36 j1260 


RESULTS AND DISCUSSIONS 
Simulation results of distribution system without AVR 


Figure 2 represents the simulation circuit of 10 bus distribution system. The simulation is carried out 
by using open modelica connection editor. A second order generator is connected as a source for the network. 
Here, 10 buses are connected in a loop structure. Bus 0 is a generator bus. The remaining buses from 1 to 10 
are load buses. To study the load flow analysis and to check stability of the network, a power fault is 
introduced in Bus 4 during the time interval of 2-3 sec. The change in voltage due to the effect of fault in 
different buses are analyzed in the simulation results. The electrical components used in this research 
simulation is taken from Open IPSL library. Table 2 represents the location of various components in 
OMEdit. 
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Figure 2. Simulation circuit of 10 bus loop distribution system without AVR 
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Table 2. Location of electrical components in OMEdit connector using [PSL 


S.No Name of the component Location of component in Open IPSL library Quantity used 
1 Generator Open IPSL-Machines-PSAT-second order synchronous Machine 1 

2 Buses Open IPSL-Electrical-Buses-Bus 11 

3 Power line Open IPSL-Electrical-Branches-PW line 10 

4 System Base Open IPSL-Electrical-System base 1 

5 Constant PQ load Open IPSL-Electrical-Loads-Constant PQ load 10 

6 Power Fault Open IPSL-Electrical-Events-PW fault 1 

7 Automatic voltage regulator Open IPSL-Electrical-Controls-PSAT-AVR-AVRTypelll 1 

8 Power system stabilizer Open IPSL-Elecrical-Controls-PSAT-PSS-PSSTypell 1 


Figure 3 represents the voltage magnitude between bus 1 and 4. Because of the fault there is a 
voltage drop in the time duration of 2-3 sec. The voltage dip in bus 1 is around 0.2pu and in bus 4 it is 0.8pu. 
After 3 sec the buses regains their actual volatges. 











Figure 3 Volatge Magnitude between Bus 1 and Bus 4 


The study of voltage between buses 9 and 4 is given in Figure 4. From the graph it is clear that the 
reduction of voltage observed in bus 9 is approximately 0.2pu. This value is lesser than the fall in voltage 
magnitude measured in bus 4. Similarly, Figure 5 represents the voltage suppression in bus 10 due to 
occurrence of fault in bus 4. The minimum voltage drop perceived is around 0.05pu which is minimum when 
compared to the voltage level observed in other buses connected in the network. Since, bus 10 is very far 
from fault introduced bus 4, the voltage suppression level during the prescribed time limit is less. 
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Figure 4. Voltage status of buses 9 and 4 


Table 3 gives a view of voltage and angle values for all the 10 buses connected in the network with 
generator and without AVR. The per unit value of voltage for the buses varies between 0.893 to 0.993. 

Also, the voltage angle varies from 101 to 108 deg. From the results, we can observe that during the 
time interval 2-3 sec all the buses experience a voltage drop before reaching its steady state value. To reduce 
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these transients an automatic volatge regualtor alongwith power system stabilizer is connected with the 
generator. It also helps to enhance the pu voltage of the buses. 





terre (s) 


Figure 5. Voltage magnitude observed in bus 10 


Table 3. Voltage magnitude and angle profile of 10 bus network without AVR 
S. No Bus number Voltage magnitude Volatge angle 


1. Bus 1 0.95771 106.183 
2 Bus 2 0.930369 104.122 
3 Bus 3 0.910809 102.576 
4 Bus 4 0.898692 101.619 
5 Bus 5 0.893786 101.299 
6 Bus 6 0.89602 101.634 
7 Bus 7 0.905495 102.615 
8 Bus 8 0.926821 104.151 
9 Bus 9 0.955941 106.2 

10 Bus 10 0.99322 108.667 


3.2. Simulation results of proposed network with AVR and PSS 

Figure 6 represents the simulation diagram of distribution network energized by a generator 
connected with AVR and PSS. A power faut is given in bus no. 4 between the time interval of 2-3 sec. A 
third order AVR and second order PSS is selected from IPSL library. 

Figure 7 gives information regarding the connection of AVR and PSS with generator. Here, V¢ of 
AVR is connected with V¢ of Generator. Connect Vto to Vr and V to V. Also, the input to PSS is taken from 
W of generator. The ouput of PSS is given as input to AVR through V, terminal. Pm and Pm0O of generator is 
interconnected. The power output of generator is given to Bus 0. 


























Figure 6. Simulation circuit of proposed network with AVR and PSS 


Fig 8 represents the change in voltage during the occurrence of fault in the network. Here, at 0.2s 
Bus 1 volatge has been reduced from 0.99 pu to 0.82 pu. So, the drop in it is nearly 0.17 pu. At 0.3s it attains 
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its original value. Bus 4 volatge has a drop of 0.63 pu. Fig 9 represents the voltage drop measured between 
buses 4 and 9. The drop in voltage of bus 9 is 0.13 pu which is minimum when compared with bus 4. Table 4 
describes the detailed summary of voltage magnitude and angle of 10 buses. The voltage value differs from 
0.978 pu-0.998 pu. Angle of the buses varies between 74.25 to 75.93 deg. 


ARTI = yi 





Figure 7. Connection diagram of AVR and PSS with generator 
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Figure 8. Voltage profile of buses 1 and 4 with AVR 
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Figure 9. Voltage drop measurement between buses 4 and 9 
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Figure 10. Voltage between buses 4 and 10 


Table 4. Outline of voltage magnitude and angle of buses with AVR 
S. No Bus Number Voltage Magnitude Volatge Angle 


1. Bus 1 0.991363 75.339 

2 Bus 2 0.985615 74.8671 
3 Bus 3 0.981533 74.5266 
4 Bus 4 0.979096 74.3208 
5 Bus 5 0.978292 74.2518 
6 Bus 6 0.979115 74.3203 
7 Bus 7 0.981572 74.5255 
8 Bus 8 0.98565 74.8661 
9 Bus 9 0.991385 75.3384 
10 Bus 10 0.998804 75.9377 


4. COMPARISON OF SIMULATION RESULTS WITH AND WITHOUT AVR 

The per unit voltage of buses 1 to 10 connected with generator is given in Figure 13. From the graph 
it is understood that the voltage drop exists in all buses between 2-3 sec. This is because of fault simulated in 
bus 4 at that time interval. So, the drop at bus 4 is more than the drop monitored in other buses. Here, the 
stable voltage is obtained at 3.2 sec. Figure 14 represents the voltage magnitude of the buses in the network 
connected with AVR. All bus volatges are interrupted at the time interval of 2-3 sec with an injected fault in 
bus no. 4. Due to this fault, the voltage drop observed in the buses like 3.5 and 6 nearer to bus 4 experience 
maximum drop. Exactly at 3 ses the fault is cleared and all the buses regains its original voltage. 

Table 5 compares various parameters between distribution network with and without AVR. The 
attributes like Voltage magnitude, Voltage drop in near and far buses, transients and stability of the system 
have been discussed here. The voltage magnitude of network without AVR varies from 0.893 to 0.993 
whereas the system with AVR varies between 0.978-0.998. Also, the drop in near bus is around 0.8pu for the 
former network and it is about 0.5pu for the later. Similarly, drop in far bus is 0.05pu in first case and 0.02pu 
in second case. Regarding the transients and stability of the network, AVR implemented network is stable 
than the other network. 
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Figure 13. Per unit Voltages spotted in all the buses without AVR 
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Figure 14. Voltage magnitude of all buses in the network with AVR 


Table 5. Comparison between power networks with and without AVR 





S. No Parameter Distribution network without AVR Distribution network with AVR 
L. Volatge magnitude 0.893-0.993pu 0.978-0.998 pu 
2. Maximum voltage drop in bus 4 0.8pu 0.6pu 
3. Maximum voltage drop in bus 5 0.8pu 0.5pu 
4. Maximum voltage drop in bus 10 0.05pu 0.02pu 
5. Transients Present Absent 
6. Stability Unstable Stable 
CONCLUSION 


Thus, from the discussion it is decided that the load flow analysis is carried out to evaluate voltage 


magnitude and angle of all load buses connected in the loop distribution network using NR method. Also, the 
voltage stability is analysed by using OMedit simulation. From the simulation results it is observed that the 
network with AVR controlled generator performs better than the generator excited system. The following 
points are observed in simulation results of AVR connected system, a) improvement in the voltage stability 
of the system, b) there is a boost up in the magnitude of the voltage, and c) drop in voltage is reduced during 
the fault injected time duration. As a future scope, the distribution network can be simulated using matlab for 
harmonic analysis to guarantee the supply of quality power to the patrons. 
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